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ABSTRACT 

We describe 2D gasdynamical models of jets that carry mass as well as energy to the 
hot gas in galaxy clusters. These flows have many attractive attributes for solving the 
galaxy cluster cooling flow problem: Why the hot gas temperature and density profiles 
resemble cooling flows but show no spectral evidence of cooling to low temperatures. Using 
an approximate model for the cluster A1795, we show that mass-carrying jets can reduce 
the overall cooling rate to or below the low values implied by X-ray spectra. Biconical 
subrelativistic jets, described by several ad hoc parameters, are assumed to be activated 
when gas flows toward or cools near a central supermassive black hole. As the jets proceed 
out from the center they entrain more and more ambient gas. The jets lose internal pressure 
by expansion and are compressed by the ambient cluster gas, becoming rather difficult to 
observe. For a wide variety of initial jet parameters and several feedback scenarios the 
global cooling can be suppressed for many Gyrs while maintaining cluster temperature 
profiles similar to those observed. The intermittancy of the feedback generates multiple 
generations of X-ray cavities similar to those observed in the Perseus Cluster and elsewhere. 
Subject headings: X-rays: galaxies - galaxies: clusters: general - X-rays: galaxies: clusters 
- galaxies: cooling flows 



1. INTRODUCTION 

Many computationally elaborate calculations of 
jet-heated cooling flows have been proposed to ex- 
plain why the hot gas in galaxy clusters cools much 
slower than the rate expected from the cluster X- 
ray luminosity, M c f w (2fj,m p /5kT V i r )L x (e.g. Cav- 
aliere et al. 2001; Reynolds, Heinz, & Begelman 
2001,2002; Bruggen & Kaiser 2002; McCarthy et 
al. 2003; Hoeft et al. 2003; Basson & Alexan- 
der 2003; Omma et al. 2004; dalla Vecchia et al. 
2004; Zanni et al. 2005). In principle the energy 
created by gas accreting at rates very much less 
than M c f onto supermassive black holes in cluster- 
centered galaxies can easily balance the radiative 
losses L x . Consequently, AGN jets have been re- 
garded as a plausible means to distribute this en- 
ergy throughout the cluster gas. Unfortunately, few 
if any jet-heated flow simulations closely resemble 
in projection the X-ray emission and radial temper- 
ature profiles observed in galaxy groups and clus- 
ters. Surprisingly often these calculations do not 
include radiative losses and therefore lack the es- 
sential physics to prove that cooling can be stopped 
by the jets. Nor do these calculations typically con- 
tinue for many Gyrs. Multi-Gyr calculations are es- 
sential to determine if the proposed type of heating 
can keep the group/cluster gas from cooling while 
also preserving time-averaged temperature and en- 
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tropy profiles commonly observed, most of which are 
very similar. In general, the long term effectiveness 
of powerful jet heating in stopping the cooling has 
not been convincingly demonstrated. As a first step 
toward a more successful solution, adopted here, jet- 
heating scenarios are sought that are consistent with 
the global observations of the hot cluster gas. Once 
this problem is solved, the next step is to determine 
if real black holes can indeed supply energy in the 
form and amount required. 

We consider here jet-dominated flows in which gas 
near the cores of cluster-centered elliptical galax- 
ies is accelerated in bipolar outflows when triggered 
by a supermassive black hole feedback mechanism. 
Both mass and energy are transported out from the 
center. The feedback appears largely as an out- 
ward flow of mass rather than energy, although the 
energetics of successful flows must be understood 
a postiori. We find that these more massive jets 
have many desirable long term attributes. We have 
discussed elsewhere how heated, buoyant gas can 
rise upstream in cooling flows, transporting mass 
and energy to distant cluster gas in a manner that 
sharply reduces the overall cooling rate and pre- 
serves the thermal profiles observed (Mathews, et 
al. 2003; Mathews, Brighenti & Boute 2004). We 
now turn our attention to the possibility that mass- 
loaded subrelativistic jets can accomplish the same 
desirable results. 

Another inspiration for the calculations described 
here are recent HI (Morganti et al. 2004; 2005), UV 
(e.g. Crenshaw et al. 1999; Kriss 2003) and X-ray 
(George et al. 1998; Risaliti et al. 2005) observa- 
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tions of AGNs showing blue-shifted absorption or 
emission lines along the line of sight. Neither the 
optical depth nor the covering factor of the outflow- 
ing gas can be accurately determined from these 
observations, but the outflowing mass flux can be 
comparable to the Eddington limit, implying that 
relatively little gas is being directly accreted by the 
central black hole. The outflow velocity is typically 
several 100 km s _1 but increases to several 1000 km 
s _1 in a few more luminous objects (Kriss 2004). At 
least half of all AGNs exhibit outflows so it is plausi- 
ble that they exist in all objects and with substantial 
covering factors. Winds from accretion disks are one 
possible explanation for the outflows (Narayan & Yi 
1994; Konigl & Kartje 1994; Blandford & Begelman 
1999; Proga 2000; Soker & Pizzolato 2005), partic- 
ularly those with higher velocities. 

The highly uncertain outflowing mass flux ob- 
served in low luminosity galaxies (M ~ 1 M Q yr -1 ), 
is far less than the mass flux that arrives at the 
central galaxies in group and cluster cooling flows, 
M ~ 10 — 300 Mq yr _1 . However, outflowing 
gas in these more massive flows is likely to be too 
hot, too rarefied and too highly ionized to produce 
blucshiftcd UV or X-ray lines and would be diffi- 
cult to observe. Indeed, currently available UV and 
X-ray absorption observations may naturally select 
the densest, coldest and most slowly moving gas in 
each outflow. Many authors (e.g. Churazov et al. 
2002; 2005; Peterson & Fabian 2005) argue that the 
mechanical outflow luminosity L mec h from massive 
black holes can greatly exceed their bolometric lu- 
minosity, and this may apply to all radiatively inef- 
ficient black holes (Hopkins et al. 2005). 

As in the Blandford-Begelman accretion model, 
we assume that the majority of gas accreted at the 
outer edges of the accretion disks ultimately flows 
away from the disk surface as a fast but nonrelativis- 
tic wind. Such disk flows are unlike the relativis- 
tic e ± radio-emitting jets driven from much smaller 
regions near the central hole (Blandford & Znajek 
1977), but the two types of outflow may coexist. In 
contrast, winds driven by radiation pressure from 
thin disks, such as those described by Proga (2000) 
are confined by strong radial, polodial fields (e.g. 
Blandford & Payne 1982), and may be directed 
mostly along the equatorial plane. However, more 
axis-oriented outflows may be possible from thicker 
disks or with different, more vertical field geome- 
tries (e.g. Everett, Konigl & Kartje 2001). Never- 
theless, fast wind-driven bipolar flows, such as we 
consider, may at their origin fill cones of significant 
solid angles. For example, Proga (2003) finds that 
most of the mass flux in disk winds at high velocities 
v > 1000 km/s, which we require in our flows, lies 
within 20 or 30 degrees from the polar axes. Such 
broad outflows are expected to entrain additional 
ambient gas as they proceed out from the central 
region. Therefore, for our exploratory calculations 
we adopt a bipolar wind geometry similar to the 



prevailing geometry in models of jet-heated cooling 
flows cited above, but allow the jets to have larger 
angular sizes at their source. Most of the mass out- 
flow in our jets arises not from the origin but by 
cntrainment of ambient gas at larger radii - such 
a model is supported by the recent observations of 
Sun, Jerius & Jones (2005) discussed below. Omma 
et al. (2004) considered the initial transient flow re- 
sulting from a bipolar outflow of this sort, but did 
not address the important question whether or not 
such jets can shut down the cooling flow for many 
Gyrs while preserving the observed thermal gradi- 
ent in the hot gas. This is our objective here. 

In the 2D gasdynamical models described here 
nonrelativistic outflows are generated by assigning 
a fixed velocity to gas that flows into a biconical 
source region (radius of a few kpc with half angle 
6j ~ 5 — 20°) at the cluster center. The acceleration 
of gas in the source bicone is activated by a feed- 
back recipe triggered as gas flows into the innermost 
zones. The 2D biconical outflows, which may for ex- 
ample represent a disk wind, proceed along the axis 
of symmetry of the computational grid. Even when 
the initial outflow has rather substantial opening 
angles (i.e. 6j <~ 20°), the flow rapidly concentrates 
within ~ 30 kpc into a much narrower jet. This 
compression occurs because the rapid pressure drop 
in the jet due to expansion causes the jet to be com- 
pressed and narrowed by the ambient gas pressure 
which decreases less rapidly with radius in the clus- 
ter gas. As the jet proceeds, it entrains additional 
ambient gas and its mass flux increases. These solu- 
tions have two excellent attributes: after many Gyrs 
the time-averaged gas temperature profile resembles 
those observed, dT/dr > in r < 0.1r vir , and very 
little gas cools below ~ T v i r . The jet itself is diffi- 
cult to observe. Finally, because of the intermittent 
nature of the feedback jet excitation, multiple gen- 
erations of large X-ray cavities are created with 2 or 
4 visible at any time, very similar to Perseus (e.g. 
Fabian et al. 2005). 

2. THE CLUSTER A1795 

We compare our gasdynamical calculations with 
the temperature and density profiles of the well- 
observed cluster Abell 1795 (Tamura et al. 201; 
Ettori et al. 2002) assumed to be at a distance 243 
Mpc. Abell 1795 is a typical rich cluster with a 
central cD galaxy and a reasonably relaxed overall 
structure (Boute & Tsai 1996). Abell 1795 has the 
usual attributes of normal cooling flows: strong cen- 
tral peak in X-ray surface brightness (e.g. Tamura 
et al. 2001), a positive temperature gradient dT/dr 
in the central regions out to ~ 500 kpc, a central 
radiative cooling time ~ 3 x 10 8 yrs that is much 
less than the cluster age (e.g. Edge et al. 1992; 
Fabian et al. 2001), optical line emission near the 
central cD (Cowie et al. 1983), an excess of blue and 
ultraviolet light possibly from massive young stars 
(Johnstone, Fabian & Nulsen 1987; Cardiel, Gor- 
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gas & Aragon-Salamanca 1998; Mittaz et a. 2001) 
and a central radio source 4C 26.42 (McNamara et 
al. 1996a, b). Chandra images near the center of 
Abell 1795 reveal an X-ray emission feature aligned 
with a remarkable optical filament (Fabian et al. 
2001). This filament and the central total mass pro- 
file, M oc r 0,6 inside 40 kpc, which is somewhat flat- 
ter than NFW (Navarro, Frenk & White 1996), may 
suggest a local deviation from hydrostatic equilib- 
rium. 

We approximate the total mass profile in A1795 
with an NFW profile with virial mass M„j r = 
10 15 Mq and concentration c = 6.57, which also 
matches the total mass found from X-ray obser- 
vations, assuming hydrostatic equilibrium. The de 
Vaucoulcurs mass profile of the central cD galaxy, 
defined by M* ~ 6 x 10 11 and R e = 8.5 kpc, has 
also been included, but this mass has little effect on 
the overall gas dynamics. 

3. COMPUTATIONAL PROCEDURE 

The two-dimensional numerical calculations de- 
scribed here are solutions of the same flow equa- 
tions described in our earlier paper on heated cool- 
ing flows (Brighcnti & Mathews 2002). These equa- 
tions explicitly include radiative cooling. The 2D 
computational grid is in spherical polar coordinates 
r and 9. Unless stated otherwise, the grid is com- 
prised of 200 logarithmically spaced radial zones ex- 
tending to 3 Mpc and having a central zone of size 
0.75 kpc. There are 60 evenly spaced angular zones 
in the range < 9 < ir. Some of our flows have been 
computed at higher resolution with 600 logarithmi- 
cally spaced radial zones extending to 3 Mpc, with a 
central zone of size 0.6 kpc, and 120 evenly-spaced 
angular zones. We find that all important results 
(temperature and density profiles, cooling rate, etc.) 
are essentially unchanged when the resolution is im- 
proved. 

When gas in a computational zone begins to cool 
by radiative losses to low temperatures, usually near 
the center of the flow, the density of cooling gas in 
the zone increases to maintain pressure equilibrium 
with gas in neighboring grid zones. This representa- 
tion of the cooling process, discretized and averaged 
on the grid scale, is unphysical since gas cooling in 
pressure equilibrium should eventually occupy a vol- 
ume much smaller than that of the grid zones. In 
order to approximately allow for this subgrid evolu- 
tion, we remove cold gas as it forms, assuming that 
its volume becomes vanishingly small. Cooling gas 
is removed by adding a mass sink term to the equa- 
tion of continuity, —q(T)p/t CO oi, where t coo i is the lo- 
cal radiative cooling time as described by Brighenti 
& Mathews (2002) and q = 2cxp(-T/T q ) 2 becomes 
large when T < T q = 5 x 10 5 K. This mass sink 
term is used to remove the unphysical clutter of 
zones containing cold gas without affecting the flow 
of hotter gas, T » T q , throughout the rest of the 
cluster. 



In some calculations unphysical cooling also oc- 
curs (generally at small radii) along the symmetry 
axis, 9 = or n, where we employ reflection bound- 
ary conditions. As gas approaches a reflecting axis, 
it is compressed and cools in a way that would not 
occur in a full 3D calculation where such reflections 
do not occur. Nevertheless, this spurious, purely 
numerical cooling near the symmetry axis is nec- 
essarily included in the computations described be- 
low. Consequently, our estimates of the cooling rate 
M may be regarded as conservative upper limits. 
When the flow velocity is entirely radial, as in the 
(unheated) cooling flow described in the next sec- 
tion, this type of boundary cooling does not occur. 

4. NORMAL COOLING FLOW IN A1795 

We begin with a simple evolutionary cooling flow 
for A1795 in which the gas is allowed to evolve from 
an initial hydrostatic model in good agreement with 
the observed density and temperature profiles. In 
this traditional spherical cooling flow all the gas 
cools at the center of the flow and there is no depen- 
dence on polar angle 9. For standard cosmological 
parameters (ft = 0.3; A = 0.7; H = 70 km s _1 
Mpc -1 ) large clusters like Abell 1795 formed rela- 
tively recently, so we consider the internal flow evo- 
lution for only 7 Gyrs. The dotted lines in Figure 
1 show the radial variation of the gas density and 
emission-weighted temperature in the cooling flow 
after 7 Gyrs. The gas density follows the observa- 
tions fairly well beyond about 50 kpc but is system- 
atically too large closer to the center. This density 
excess is typical for pure cooling flows, as discussed 
by Mathews & Brighcnti (2003). The depression 
of the observed (azimuthally averaged) density rel- 
ative to this flow within 50 kpc may be due in part 
to unresolved X-ray cavities. 

The cooling rate M(t) for this cooling flow, shown 
as a dotted line in Figure 2, increases with time, 
approaching ~ 400 Mq yr _1 after 7 Gyrs. This 
cooling rate is comparable to the cooling rate for 
A1795 estimated by Allen et al. (2000) from depro- 
jected ROSAT images, M ~ 500 M Q yr" 1 . How- 
ever, XMM RGS spectra show no evidence of gas 
with temperatures less than ~ 2 keV (Tamura et 
al. 2001), indicating a much smaller cooling rate, 
M < 150 Mq yr _1 . This upper limit is consistent 
with M < 100 Mq yr _1 estimated from Chandra 
observations (Ettori et al. 2002). 

5. JET OUTFLOWS 

Gas flows including the effects of jet momentum 
are solved in several stages. Each calculation begins 
with a static cluster atmosphere based on the ob- 
served temperature and density profiles in A1795. 
During time < t < 1 Gyr the configuration is 
allowed to evolve (without jets) toward a pure cool- 
ing flow, and later during 1 < t < 7 Gyrs the jet 
momentum is activated according to a feedback cri- 
terion. 
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We adopt a simple computational procedure to 
trigger jet outflow in which all gas in a biconical 
jet source region near the center of the flow is set 
into outward motion at velocity Uj as long as some 
feedback criterion is satisfied. The geometrical pa- 
rameters that define the jet source region are the 
radius rj and the half opening angle 9j of the jet. 
We consider three feedback criteria to activate the 
jet: 

A: The gas velocity in the jet source re- 
gion is set to Uj only when the gas cool- 
ing rate M is non-zero inside a radius of 
1 kpc. 

B: The gas velocity in the jet source re- 
gion is set to uj only when the net mass 
flow across a sphere of radius 1 kpc is 
negative (mass inflow). 

C: Continuous jet flow Uj in the jet 
source region at all times. 

During times when the jet outflow is not active, gas 
flows through the source region in accordance with 
the usual gasdynamical equations. Each flow calcu- 
lation is uniquely designated by mN(X,rj,Oj,Uj) 
where N is a number assigned to each computed 
flow, X = A, B or C is the feedback criterion, rj is 
the jet radius in kpc at the source, 9j is the source 
jet half-angle in degrees and Uj is the jet source ve- 
locity in units of 10 3 km s" 1 . Flows at the higher 
spatial resolution are designated with upper case 
MN{X,rj,6j,Uj). 

5.1. A Representative Flow with Jet Momentum 

Among the many models with satisfactory 
or excellent results, we select ml(A, 5, 10, 10) 
as representative and discuss it in more de- 
tail. The azimuthally averaged gas density and 
emission-weighted temperature profiles for the 
m\(A, 5, 10, 10) flow are shown at three times in 
Figure 1. The global cooling rate M(i) for this 
flow, plotted in Figure 2, is small, indeed its time- 
averaged value (M(t)) 20 M© yr _1 is well be- 
low the constraints imposed by XMM and Chan- 
dra observations. Even more remarkable, both the 
density and temperature profiles shown in Figure 
1 retain their cooling flow appearance at time 7 
Gyrs. Within about 50 kpc n(r) and T(r) for the 
m(A, 5, 10, 10) flow lie between the pure cooling flow 
and the observations. 

Figure 3 shows in more detail the 2D density 
structure at several times for the high resolution ver- 
sion of this flow, Ml (A, 5, 10, 10). The four panels 
in Figure 3 show illustrate the growth of successive 
generations of buoyant X-ray cavities. Cavities are 
associated with jet intermittancy. The evolution of 
the buoyant cavities away from the (horizontal) jet 
axis may be an artifact of the 2D reflection bound- 
ary conditions along this axis, but it is comparable 



to a similar deviation found in the 3D models of 
Bruggen et al. (2002) . Nevertheless, the creation of 
multiple pairs of bubble cavities, similar to those in 
Perseus (Fabian et al. 2003), is an encouraging fea- 
ture of these jet flows. Finally, we stress that these 
jets carry mass as well as energy to large distances 
from the central AGN and in this respect they differ 
from many previous calculations in which the jets 
carried little or no mass. 

5.2. Additional Jet Momentum Flows 

Table 1 summarizes some of the jet-heated flows 
we have calculated. For consistency, all results in 
Table 1 refer to flows computed at lower resolution, 
but are not significantly changed when the grid res- 
olution is refined. In addition to listing the param- 
eters that define each flow, Table 1 gives several ad- 
ditional global results after 6 Gyrs of jet feedback: 
the total energy Ekin supplied by the jet source, 
the time-averaged mechanical luminosity generated 
in the jet source region L mec /i, the total mass that 
cooled M coo i after t = 7 Gyrs and the average cool- 
ing rate (M), and the X-ray luminosity L x at this 
same time. For many jet parameter combinations, 
the mean cooling rate (M) is less than or compara- 
ble to the Chandra value M w 100 M & yr" 1 . The 
mechanical energy for successful flows ranges from 
~ 0.02L X to <~ L x . We stress that the condition 
Lmech ^ L x that is usually required to keep jet- 
heated flows from cooling does not necessarily ap- 
ply to our flows in which mass as well as energy is 
transported outward. Consequently, successful jet- 
advecting flows are possible even when L mec h -C L x . 
Jet mass transfer appears to be an efficient and ro- 
bust way to recirculate gas and energy outward with 
little radiative cooling while retaining the cooling 
flow appearance as observed in the density and tem- 
perature profiles. 

Regarding this latter important point, Figure 4 
shows the gas density and temperature at time t = 7 
Gyrs for a sample of successful and unsuccessful 
flows listed in Table 1. The cooling rates for these 
flows (sampled each 0.5 Gyr) are shown in Figure 
5. In both Figures 1 and 4, the computed gas den- 
sity within about 20-30 kpc from the center exceeds 
that observed in A1795, but this region of A1795 
contains a cool, transient optical filament about 60 
kpc in diameter and may be experiencing a local 
deviation from hydrostatic equilibrium. For many 
of the unsuccessful models with (M) > 100 Mq 
yr _1 , most of the cooling occurs during one or two 
episodes. These flows would be regarded as suc- 
cessful in matching the X-ray data if they were com- 
puted for only ~ 2 — 3 Gyrs when M(t) is acceptably 
small, but star formation, which is likely to occur 
during times of enhanced cooling, would be inconsis- 
tent with optical colors of cluster-centered galaxies 
(e.g. McNamara 1997). Enhanced blue light from 
young stars typically persists for 2 - 3 Gyrs after a 
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starburst. 

The mass M coo i in column (9) of Table 1 rep- 
resents all the cooled gas that has been removed 
from the grid by the term —q(T)p/t coo i in the equa- 
tion of continuity. Almost all of this gas cools in 
the central regions. However, M coo i is not small 
and is often much larger than the mass of central 
black holes (or stars) in cluster-centered galaxies. 
However, there are reasons to believe that we have 
overestimated the global cooling rate in our mod- 
els. For example the sink term —q(T)p/t coo i in the 
continuity equation not only removes the cooled gas, 
but may also encourage local pressure gradients that 
stimulate additional cooling. If gas in a particular 
grid zone begins to cool while gas is being removed 
by the sink term, hot gas from adjacent zones, that 
would not otherwise cool, may be stimulated to flow 
toward the cooling zone, possibly raising the local 
cooling rate unrealistically. As discussed above, we 
also expect the computed cooling rate to be spuri- 
ously enhanced by cooling near the symmetry axis 
where reflecting boundary conditions must be em- 
ployed. Evidence for computational overcooling is 
provided in the final two models listed in Table 
1, ml3(B,5, 10,10) and ml4(B,5,10,5), in which the 
sink term is set to zero, q(T) = 0. In flows with 
q = 0, M coo i in Table 1 represents the mass of cooled 
gas (T « T q = 5 x 10 5 K) that remains in the 
grid and goes into approximate free fall if it is not 
at the center. Values of M coo i for these flows arc 
very small, suggesting that nonzero q(T) does in- 
deed artificially increase the cooling rate. But real 
cooling at some level can and does occur. Signifi- 
cant centrally-located cold gas and star formation 
are observed in many massive clusters (e.g. Edge 
2001). Recent CO observations of A1795 by Salome 
& Combes (2004) have detected ~ 10 11 M of cold 
gas which is entirely consistent with our flow calcu- 
lations based on an approximate model for A1795. 

In marked contrast to our earlier models in which 
we explored a wide variety of heated cooling flows 
(Brighenti & Mathews 2002; 2003), very few of the 
jet momentum flows described here have tempera- 
ture and density profiles that strongly deviate from 
the observations of A1795. The relative success of 
each jet momentum flow, expressed in column (11) 
of Table 1, is based largely on the magnitude of 
the mean cooling rate (M). The results in Ta- 
ble 1 and Figures 2 and 5 imply limits on the jet 
source parameters X, rj, 6j, and uj corresponding 
to (M) < 100 Mq yr" 1 as observed with Chan- 
dra. In particular, radiative cooling is effectively 
shut down for source regions with radii rj > 3 kpc, 
half angles 9j > 10° and jet velocities Uj > 5000 km 
s- 1 . 

5.3. Nature of the Jet Flow 

The simple jets we employ are based on the 
plausible notion that strong non-relativistic winds 
flow from accretion disks around supermassive black 



holes and that these winds return most of the mass 
inflow received from centrally cooling gas. We 
implement this simple idea by insisting that the 
outflow velocity remains constant and rather large 
throughout the biconical jet source region when- 
ever a feedback criterion is satisfied. Although this 
model for the jet source is admittedly ad hoc, as 
the jets move further out they appear to develop 
more universal properties. In this section we briefly 
review the evolution of jets far beyond the source 
region. 

The physical nature of our jets is most clearly 
defined when the jet activity is continuous in 
time (models C) such the high resolution model 
M11(C, 5,20,10) in which the overall flow ap- 
proaches a quasi-steady state. The global velocity 
field and density contours for the M11(C, 5, 20, 10) 
flow after 7 Gyrs are shown in Figure 6. The mag- 
nitude of the radial gas flow at time 7 Gyrs in 
three angular directions are shown in Figure 7. It is 
apparent from this Figure that the cooling inflow, 
which fills most of the cluster volume, flows toward 
the jet, becomes entrained and is carried outward. 
The continuous jet creates an extended two-sided 
channel of low density gas. The difficulty of ob- 
serving such a channel can be seen in the X-ray 
surface brightness images of the (high resolution) 
M1(A, 5, 10, 10) and MU(C, 5, 20, 10) flows shown 
in Figure 8, viewed perpendicular to the jet. From 
this viewing direction the jet cavitation in the con- 
tinuous M11(C, 5, 20, 10) flow produces at most a 
10% reduction in the X-ray surface brightness along 
the jet axis. The jet cavitation produced by inter- 
mittent flows such as M1(A, 5, 10, 10) is less pro- 
nounced and more difficult to detect. Although our 
2D jets are constrained to flow in a axisymmetric 
fashion along the 9 = and it axes, in a full 3D 
simulation the jets may not follow such a perfectly 
linear pathway through the cluster gas. 

The symmetric jet cavitations visible in Figure 8 
are remarkably similar to the Chandra observations 
of the double-jetted radio galaxy NGC 1265 by Sun, 
Jerius & Jones (2005). The X-ray contours in their 
Figure la show symmetric indentations along the 
radio jet axis just as in our Figure 8, suggesting that 
ambient hot gas in NGC 1265 is being entrained 
and swept along with the jet. This observation also 
supports the biconical outflow geometry that we as- 
sume here rather than disk winds in the equatorial 
plane (e.g. Proga 2000). 

It is apparent from Figure 8 that our jets are 
poorly resolved, often occupying only a few angu- 
lar zones even at the higher grid resolution. This 
poor resolution restricts somewhat our analysis of 
the transverse jet profiles. Nevertheless, the behav- 
ior of the gas outside of the small jet regions - which 
determines the overall X-ray properties of the clus- 
ter - does not appear to be strongly affected by the 
level of numerical resolution in the jets. 

To understand better the details of the jet- 



6 



atmosphere interaction, we have examined the high 
resolution jet in the M11(C, 5, 20, 10) flow. One 
natural attribute of the biconical jet source region 
is that the original gas within the source is expelled 
very rapidly after the jet turns on. Consequently, 
during most of the active phase of this region, most 
of the outgoing gas is supplied by recent advection 
near the outer boundary of the bicone. As a re- 
sult, most of the mass outflowing from the jet at rj 
occurs near 6j so the initial jet density profile per- 
pendicular to the jet axis is hollow with a strong 
central minimum. This type of hollow jet structure 
may in fact be physically appropriate if entrainment 
of ambient gas occurs near the source region. 

An important feature of these jet solutions is 
the narrowing of the jets as they moves outward. 
Although the half angle of biconical jets 9j can 
be rather large, these initially conical jets become 
nearly cylindrical as they move outward. This jet 
focusing occurs because the pressure in the jet is 
rapidly lowered by expansion and tends to decrease 
with radius faster than the pressure in the ambi- 
ent gas. As a result the jets are compressed and 
collimated by the ambient gas pressure. 

A number of (poorly resolved) internal shocks ap- 
pear in the jet flow that we do not describe in detail 
here. The main effect of these shocks, even inside 
the jet source region, is to raise the temperature 
and entropy within the central jet to rather high 
values. Gas in the jet source and just beyond flows 
approximately at the sound speed within the jet. 
As the jet moves further out, the mass outflow in 
the jet increases due to entrainment. The amount 
of entrainment is approximately independent of the 
spatial resolution. The top panel in Figure 9 shows 
the variation of the (outward) mass flux transverse 
to the jet axis in the continuous M11(C, 5, 20, 10) 
jet at various cluster radii. The quantity plotted is 
dM/dO = pv r 2irr 2 sin(0) M Q yr" 1 where v r (6) is 
the local jet velocity. It is clear from Figure 9 the 
angular width of the jet narrows from its initial half 
angle 6j — 20° as it moves out and remains hollow 
to rather large distances in the cluster. 

Of particular interest is the radial increase of 
the integrated mass flow in the jet, M(9) = 

Jq (dM /d9)d9, plotted in the lower panel of Figure 
9. The total mass outflow in the jet and its ap- 
proximate local angular width 9 max (r), can be es- 
timated from the maximum M = M(0 max ) in each 
curve. The decline in M{9) for 9 > 9 max is due to 
the negative contribution to the integrated flow be- 
yond 9 m ax caused by slowly inflowing gas adjacent 
to the jet. The mass flux increases from 65 M /yr 
at 5 kpc to 116 M Q /yr at 20 kpc and reaches 170 
M /yr at 200 kpc. This increase shows that mass 
entrainment is a key feature of the success of these 
simulations. For model M11(C, 5, 20, 10) in which 
the jet is continuously active, M{9 = 90°) becomes 
essentially zero, indicating that the jet is returning 



mass to large radii at the same rate that it arrives at 
the center in the cooling inflow outside the jet. This 
mass conservation is approximately true for all the 
other flows, for example M(90°) « -12 ± 4 M Q /yr 
for the flows plotted in Figure 4. 

We find that the decelerating jets penetrate to 
large distances in the cluster gas, well beyond sev- 
eral hundred kpc. However, the maximum distance 
to which the jet outflow continues is observed to 
increase with the refinement of the computational 
grid and this must be explored in future calcula- 
tions. In fully 3D versions of the MU(C, 5, 20, 10) 
jet we expect that the lateral motions of jet due to 
shear instabilities may cause the jet to dissipate its 
energy at a somewhat smaller radius. 

6. CONCLUSIONS 

Subrclativistic jet flows that entrain ambient gas 
may be the essential key for solving the cluster cool- 
ing flow problem: Why the hot gas temperature and 
density profiles resemble cooling flows but show no 
spectral evidence for cooling below ~ T vir /3 at rates 
expected from the luminosity L x . Many scenarios 
have been previously considered in which the gas is 
heated by a variety of mechanisms including jets. In 
most previous jet calculations it has been assumed 
that the jets are primarily sources of energy that 
reheat the cluster gas with no significant outward 
mass transport. Overall, these simulations have not 
been successful in reducing the cooling rate while 
maintaining the observed temperature and density 
profiles for many Gyrs. Because of the large num- 
ber of cooling core clusters observed, cooling must 
be sharply reduced or arrested for many Gyrs. 

The mass-carrying jets considered here are a vari- 
ant of the circulation flows we have discussed in 
which buoyant bubbles provide an outward mass 
transport while most of the X-ray emission comes 
from a normal cooling interbubble inflow (Mathews, 
et al. 2003; Mathews, Brighenti & Boute 2004). 

Our computational results for mass-carrying jets 
are robust in that wc find satisfactory multi-Gyr 
solutions for a significant range of parameters de- 
scribing the initial central jet outflow. The impor- 
tant global features of our flows are also insensi- 
tive to computational resolution. The jet outflow 
is stimulated by cooling or inflowing gas near the 
central supermassive black hole, but the success of 
our longterm solutions is not strongly dependent on 
the specifics of this feedback mechanism. The nat- 
ural intcrmittancy of the feedback generates mul- 
tiple generations of X-ray cavities similar to those 
observed in the Perseus Cluster and elsewhere. Nev- 
ertheless, the physics at the source of outflow is 
poorly understood and observational support is lim- 
ited. More detailed models of mass-carrying jets will 
be necessary before they can be fully accepted. 

One possible objection to the jet driven mass cir- 
culation described here is that SNIa-enriched gas 
that enters the source cone (within the central E 
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galaxy) is transported out only along the jet axis, 
unlike the observed iron abundance pattern which 
is spherically symmetric around the central galaxy 
(de Grandi et al. 2004). Either the radio jet pro- 
cesses (Gower, et al. 1982) as in the cluster observed 
by Gitti et al. (2005) or the jet axis direction was 
altered by black hole mergers at early times when 



most of the SNIa iron was produced. 

Studies of the evolution of hot gas in elliptical 
galaxies at UC Santa Cruz are supported by NASA 
grants NAG 5-8409 & ATP02-0122-0079 and NSF 
grant AST-0098351 for which we are very grateful. 
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TABLE 1 
GASDYNAMICAL MODELS 



model 


feed- 


r i 


°j 


Uj 

(10 3 
km s ~~ ) 


k xti 




L x 


M coo l 


<M> 


comment 3 




back 


(kpc) 


(°) 


(10 62 
ere) 


(10 45 
ere: s — ) 


(10 45 
ere s — ) 


(10 11 

Mr,) 


(M Q yr- 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


ml(A,5,10,10) 


A 


5 


10 


10 


4.90 


2.59 


2.19 


1.20 


20.0 


OK 


m2(A,5,20,10) 


A 


5 


20 


10 


2.43 


1.28 


2.63 


0.99 


16.56 


OK 


m3(A,5,20,5) 


A 


5 


20 


5 


1.45 


0.77 


3.23 


3.72 


62.1 


marginal 


m4(B,5,10,10) 


B 


5 


10 


10 


6.71 


3.54 


2.19 


1.50 


25.1 


OK 


m5(B,5,10,5) 


B 


5 


10 


5 


1.16 


0.61 


2.87 


9.14 


152.4 


fails 


m6(B,5,20,10) 


B 


5 


20 


10 


5.54 


2.93 


2.17 


0.016 


0.27 


OK 


m7(B,5,20,5) 


B 


5 


20 


5 


2.10 


1.11 


2.65 


0.47 


7.88 


OK 


m8(B,5,20,l) 


B 


5 


20 


1 


0.11 


0.054 


2.31 


21.5 


359 


fails 


m9(A,3,20,10) 


A 


3 


20 


10 


2.80 


1.48 


2.62 


3.81 


63.5 


marginal 


mlO(A,3,20,5) 


A 


3 


20 


5 


1.67 


0.88 


2.83 


6.22 


104 


fails 


mll(C,5,20,10) 


C 


5 


20 


10 


0.087 


0.046 


2.18 


0.016 


0.27 


OK 


ml2(C,5,20,5) 


C 


5 


20 


5 


0.40 


0.21 


2.65 


0.47 


7.88 


OK 


ml3(B,5,10,10) b 


B 


5 


10 


10 


7.65 


4.09 


2.11 


~ 0.02 


~ 


OK 


ml4(B,5,10,5) b 


B 


5 


10 


5 


7.36 


3.88 


1.94 


~ 


~ 


OK 



a Evaluation of the relative success of each computed flow after t — 7 Gyrs is based primarily on requiring (M) < 50 Mq yr 
b In these models q{T) = 0. 
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Log r (kpc) 



Fig. 1. — The observed hot gas density and temperature in A1795 arc shown with filled triangles (XMM observations 
from Tamura et al. 2001) and open triangles (Chandra observations from Ettori et al. 2002). The dotted line shows the 
quasi-steady pure cooling flow at time t = 7 Gyrs. The other lines show the computed density and temperature profiles for flow 
ml(A, 5, 10, 10) at three times: 2 Gyrs (long-dashed lines), 4 Gyrs (short-dashed Hnes), and 6 Gyrs (solid lines). 
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Fig. 2.— 



The total cooling rate M(t) for the cooling flow solution (dotted line) and for flow ml(A, 5, 10, 10) (solid line). 
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Fig. 3. — The density distribution for flow calculation M1(A, 5, 10, 10) at four times. 
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Fig. 4. — Density and temperature profiles at time t = 7 Gyr for six computed flows: m3(A,5,20,5) (solid line), m9(A,3,20,10) 
(long-dashed line), m7(B,5,20,5) (short-dashed line), mlO(A,3,20,5) (short dashed-dotted line), ml2(C,5,20,5) (long dashed- 
dotted line), and m7(B,5,20,5) (dotted line). 
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Fig. 5. — Total cooling rates M(t) for six computed flows: m3(A,5,20,5) (solid line), m9(A,3,20,10) (long- dashed line), 
m7(B,5,20,5) (short-dashed line), mlO(A,3,20,5) (short dashed-dotted line), ml2(C,5,20,5) (long dashed-dotted line), and 
m7(B, 5,20,5) (dotted line). 
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Fig. 6. — Velocity flow superimposed on gas density contours for the central region of computed flow M11(C, 5, 20, 10) shown 
in cylindrical coordinates. The jet becomes narrower as it moves out. 
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Fig. 7. — A snapshot of the radial gas velocity profiles in model M11(C, 5,20, 10) at time t = 7 Gyrs along three angular 
directions: the axial jet outflow at $ = (solid line), in which the plotted velocities have been reduced by 10~ 2 , 9 = n/4 
(dashed line), and the equatorial flow at 6 = tt/2 (dotted line). 
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Fig. 8. — Contours of the bolometric X-ray surface brightness for flows M1(A,5,10,10) (top panel) and M11(C,5,20,10) 
(bottom panel) at time t = 7 Gyrs. 
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Fig. 9. — Angular variation of the jet flow dM/d8 (top panel) and the integrated jet flow M(0) to (half) angle 9 (bottom 
panel) computed for the continuous jet model M11(C,5,20,10). Both profiles are shown for seven cluster radii: 5 kpc (solid 
line) 10 kpc (dotted line) 20 kpc (short dashed line) 30 kpc (long dashed line) 50 kpc (short dashed-dotted line) 100 kpc (long 
dashed-dotted line) and 200 kpc (short and long dashed line). 
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